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We perform ab initio no-core shell-model calculations for A — IS and 19 nuclei in a AhUl, or 
iVmax = 4, model space by using the effective JISP16 and chiral N3LO nucleon-nucleon potentials 
and transform the many-body effective Hamiltonians into the OhLl model space to construct the 
^-body effective Hamiltonians in the sd-shelL We separate the ^-body effective Hamiltonians with 
^ = 18 and A = 19 into inert core, one-, and two-body components. Then, we use these core, 
one-, and two-body components to perform standard shell-model calculations for the ^ = 18 and 
^ = 19 systems with valence nucleons restricted to the sd shell. Finally, we compare the standard 
shell-model results in the OhD model space with the exact no-core shell model results in the AhLl 
model space for the A = IS and A = 19 systems and find good agreement. 

PACS numbers: 21.60.De, 21.60.Cs, 21.30.Fe, 27.20.+n 


I. INTRODUCTION 

In recent years remarkable progress in ab initio micro¬ 
scopic nuclear structure studies has been made in calcu¬ 
lating nuclear properties, e.g., low-lying spectra, transi¬ 
tion strengths, etc., in light nuclei. Large-basis ab initio 
no-core shell-model (NCSM) calculations, which provide 
the foundation for this investigation, have been successful 
in reproducing the low-lying spectra and other properties 
of nuclei with < 16 [ll-llBj. 

In NCSM calculations all nucleons in the nucleus are 
active and treated equivalently in the chosen model 
space. When we increase the model space to obtain more 
precise results, we encounter the problem that the size 
of the calculations can easily exceed currently available 
computational resources. This is especially true as one 
proceeds towards the upper end of p-shell nuclei and be¬ 
yond. The problem may be cast as a challenge to re¬ 
produce the many-body correlations present in the large 
space in a tractable, smaller model space. Success in this 
endeavor will open up the prospects for ab initio solutions 
for a wider range of nuclei than are currently accessible. 

The NCSM has proven to be an ab initio microscopic 
nuclear structure approach that has been able to repro¬ 
duce experimental results and to make reliable predic¬ 
tions for nuclei with H < 16. These successes motivate 
us to develop approaches for heavier-mass nuclei. In one 
approach, a small model space effective interaction has 
been constructed by modifying the one-body piece of the 
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effective two-body Hamiltonian and employing a unitary 
transformation in order to account for many-bo dycor re¬ 
lations for the A-body system in a large space [^. In 
another approach [U, the effective two- and three-body 
Hamiltonians for p-shell nuclei have been constructed by 
performing 12hUl ab initio [i.e., A^max = 12 harmonic 
oscillator (HO) quanta above the minimum required] 
NCSM calculations for H = 6 and A = 7 systems and 
explicitly projecting the many-body Hamiltonians onto 
the OhUt space. These H-dependent effective Hamiltoni¬ 
ans can be separated into core, one-body, and two-body 
(and three-body) components, all of which are also A- 
dependent 0. 

Recently, two more ab initio methods for valence nu¬ 
cleon effective interactions have been introduced with 
the same goals; one is based on the in-medium similar¬ 
ity renormalization group approach a nd the other is 
based on the coupled-cluster method @ . 

In this work, following the original idea of Refs. 
PHI , we derive two-body effective interactions for the 
sd shell by using AhUl NCSM wave functions at the two- 
body cluster level, which contain all the many-body cor¬ 
relations of the AhUl no-core model space. The goal of 
this work is to demonstrate feasibility of this approach 
in the sd shell, where we do not require calculations at 
the limit of currently accessible computers. Such a major 
extension will be addressed in a future effort. 

At the first step, we construct a “primary” effective 
Hamiltonian following the Okubo-Lee-Suzuki (OLS) uni¬ 
tary transformation method [234^. We indicate this 
first step schematically by the progression shown with 
the two large squares in the lower section of Fig. [TJ We 
elect to perform this first step at the two-body cluster 
level for ^^F in the AfiUl model space (the “P-space”) fol¬ 
lowing the NCSM prescription PPP. For our initial 
interactions we select the JISP16 and chiral N3L0 
[HI potentials. Our formalism may be directly adapted 



2 


P 

Q 



the and systems with the primary effective 

interaction in the ATiVt basis space. The results of these 
calculations produce, respectively, the core and one-body 
components included in the secondary effective Hamito- 
nian. 

At the fourth step, we subtract the core and one-body 
terms from the secondary effective Hamiltonians of step 2 
to obtain the effective valence interaction two-body ma¬ 
trix elements (TBMEs) in the sd-shell space. 

Following the completion of these four steps, we then 
use the effective valence interaction matrix elements 
along with the extracted single-particle energies (for both 
the proton and the neutron) for standard shell-model 
(SSM) calculations in the sd-shell space. 


FIG. 1: (Color online) Flow of renormalizations adopted to 
derive an effective interaction for valence nucleons. The OLS 
procedure is first applied to derive a NCSM effective inter¬ 
action for the full A-nucleon system resulting in the “pri¬ 
mary” effective Hamiltonian PHgffP for the chosen no-core 
basis space (the “P-space”) indicated on the large square on 
the right of the figure in its upper left corner. The many- 
body truncation is indicated by iVmax, the total number of 
HO quanta above the minimum for that system. The OLS 
procedure is applied again by using the results of the NCSM 
calculation to derive the “secondary” effective Hamiltonian 
P'PgffP' for the chosen valence space (the P^space with the 
smaller many-body cutoff iV^ax) indicated on the square in 
the upper right of the figure. 


to include the three-nucleon force (3NF) but the compu¬ 
tational effort increases dramatically. Thus, we do not 
include the 3NF in this initial work. 

For the second step, we begin by performing a NCSM 
calculation for ^^F with the primary effective Hamilto¬ 
nian in the ATiQ. model space to generate the low-lying 
eigenvalues and eigenvectors needed for a second OLS 
transformation as indicated by the flow to the upper right 
in Fig. [TJ These ^^F eigenvectors are dominated by con¬ 
figurations with an system in the lowest available HO 
orbits and two nucleons in the sd shell. All additional 
many-body correlations are also present. With these ^^F 
eigenvectors and eigenvalues we then solve for the “sec¬ 
ondary” effective Hamiltonian, again using an OLS trans¬ 
formation, that acts only in the ^ space of ^^F 

but produces the same low-lying eigenvalues. Here we 
are following the scheme initially introduced in Ref. ini- 
The matrix elements of this secondary effective Hamilto¬ 
nian have the property that all configurations are defined 
with two nucleons in the sd-space and an ^^O subsys¬ 
tem restricted to the lowest available HO single-particle 
states. This second step therefore produces a secondary 
effective Hamiltonian that is equivalent to what we would 
call the 18-body cluster Hamiltonian in the NCSM acting 
in the = 0 space. 

At the third step, we carry out NCSM calculations for 


For any system with A > 18, we can obtain its 18-body 
cluster Hamiltonian by repeating the entire procedure 
utilizing the primary effective Hamiltonian for that value 
of A. The second and subsequent steps remain the same. 
That is, we perform NCSM calculations with the primary 
(A-dependent) effective NN potential for ^^O, ^^O, ^^F, 
and ^^F in order to obtain the (A-dependent) core energy, 
single-particle energies, and TBMEs, which can then be 
used in a SSM calculation for that value of A. We provide 
details for applications to A > 18 systems below using 
^^E as an example. 

We employ the Coulomb interaction between the pro¬ 
tons in the NCSM calculations which gives rise to the ma¬ 
jor shift between the derived neutron and proton single¬ 
particle energies. Exploration of full charge-dependence 
in the derived two-body valence interactions will be ad¬ 
dressed in a future effort. In particular, our current 
A = 18 and 19 applications will have at most one va¬ 
lence proton so we do not require a residual Coulomb 
interaction between valence protons in this work. 

Eor the chiral N3LO we retain full charge dependence 
in the first step that is, when deriving the primary effec¬ 
tive Hamiltonian. Thus, the A-body, core, and valence 
system calculations are performed with full charge depen¬ 
dence retained. Since we currently solve only for ^^E in 
step 2, we derive only the isospin-dependent but charge- 
independent secondary effective Hamiltonian. To retain 
full charge dependence in the secondary effective Hamil¬ 
tonian, which would constitute predictions beyond con¬ 
ventional phenomenological interactions, would require 
additional and ^^Ne calculations that are intended 
in future efforts. 

One may straightforwardly generalize these steps out¬ 
lined above to solve for effective three-body valence in¬ 
teractions suitable for SSM calculations. Earlier efforts 
using an alternative implementation of step 3 [111 showed 
that effective three-body valence interactions lead to sig¬ 
nificant improvements over effective two-body valence in¬ 
teractions. 
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II. THEORETICAL DESCRIPTION 
A. No-core shell model and effective interaction 


The NCSM calculations start with the intrinsic Hamil¬ 
tonian of the A-nucleon system, omitting any 3NF in the 
present effort: 


E 

i<j=l 


{Pi-Pjf 


2Am 


E^* 

l<j = l 


NN 


( 1 ) 


where m is the nucleon mass, is the bare NN in- 

teraction, T^ei is the relative kinetic energy and is 

the total two-body interaction. We will add the Coulomb 
interaction between the protons at a later stage since we 
treat it as a perturbative correction to the derived pri¬ 
mary effective Hamiltonian. In order to facilitate con¬ 
vergence, we modify Eq. (HD by adding (and later sub¬ 
tracting) the center-of-mass HO Hamiltonian which in¬ 
troduces a dependence on the HO energy, TiVt^ and this 
dependence is denoted by in what follows. In ad¬ 
dition, we introduce a < A to define a new a-, A-, and 
f^-dependent Hamiltonian: 


Ha = Y. 


Pi I o2 ^2 
2m 2 


a 

+ ^ Vij{n,A), ( 2 ) 

i<j=l 


Hamiltonian H 2 in the T = 0 channel different from the 
deuteron Hamiltonian. In order to preserve Galilean in¬ 
variance in the primary effective Hamiltonian, we obtain 
the solutions to Eq. (|4|) in the relative HO basis where the 
the center-of-mass component of the first term in Eq. m 
plays no role. 

We now introduce our representation of the unitary 
transformation needed to construct the primary effective 
Hamiltonian PIpQf^P := H 2 in the P space (signified 
by a superscript “P”) of the first step. The P-space 
effective interactions have A dependence, Cl dependence 
and Amax dependence all implied by the superscript P. 
We define A^max as the maximum number of HO quanta in 
the many-body HO basis space (the NCSM basis space) 
above the minimum for the A-nucleon nucleus. We select 
A^max = 4 in the present work. The resulting finite P 
space, of dimension dp, for the first step is indicated on 
the left-hand side of Eig. [TJ The diagonalization of the 
Hamiltonian H 2 in the relative HO basis provides the 
unitary transformation U 2 such that 

= (5) 

where P2;diag is the diagonal matrix containing the eigen¬ 
values p 2 ;/e- 


Ho. 


2;diag 


p 2 ;l 0 • • • 0 

0 p 2;2 • • • 0 

0 0 0 p2;max 


(6) 


where a = A corresponds to the full Hamiltonian of 
Eq. ^ with the center-of-mass HO Hamiltonian added 
and Vij{Cl, A) is the modified bare NN interaction which 
we define independent of the parameter a but including 
dependence on A: 

Vij {n, A) = - ^(rl - fjf. (3) 


where the subscript “max” signifies the dimension of the 
a = 2 space sufficient to guarantee convergence of the 
dp low-lying eigenvalues and eigenvectors. We typically 
employ max = 200 to 450 for a realistic NN interac¬ 
tion governed by the need to converge the results for the 
chosen interaction at the selected value of TiQ. 

By introducing the model space P, one builds the ma¬ 
trix Pl^diag ~ ^^2;diag^- 


The exact solution of Eq. (HD for a subset of its eigen- 
solutions in a finite matrix diagonalization requires the 
derivation of an A-b ody effective interaction for suffi¬ 
ciently heavy nuclei [I 3 , but such a derivation is not 
currently possible for A > 5 with realistic interactions. 

Here, we adopt the two-body cluster approximation 
(a = 2) for the effective interaction [mil. This allows 
us to solve the eigenvalue problem for a sufficiently large 
basis space that we achieve convergence of a suitable set 
of low-lying eigenvalues and eigenvectors needed to con¬ 
struct the primary effective Hamiltonian. In the a = 2 
approximation, the Hamiltonian (ED becomes 


ff2 = E 


2m + 2 * 


+ Vi2{^,A). 


(4) 


/ E^ 2;1 0 ... 0 \ 

", ::: ", ■ m 

\ 0 0 0 E2,dpl 

The unitary transformation matrix Ua in which a = 2 
refers to two-body cluster approximation and can be split 
into four blocks corresponding to the blocks within the 
large squares of Eig. [TJ 

/ TjP tjPQ \ 

“ ( ur vs) ■ <«> 

where the matrix is the dp x dp square matrix corre¬ 
sponding to the P space. One constructs the matrix 
from the Ua matrix by taking dp rows and columns of the 
eigenvectors corresponding to the chosen dp eigenvalues: 


Eor deriving an effective three-nucleon interaction one 
would take a = 3. Note that the A-dependence enters the 
Hamiltonian H 2 through the second term in Eq. (j3j) . Eor 
example, this A-dependence makes the two-body cluster 


U2^ 


bi,i 

^1,2 • 

• • h,dp \ 


0 - 

to 

^2,2 • • 

• ^ 2 ,dp 

(9) 

V bdp,i 

^dp,2 

• • ^dp,dp / 
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The primary effective Hamiltonian H 2 ^ signified by the 
box labeled in Fig. [TJ can then be calculated 

by using the following formula: 

TjP 

TjP _ ^2 TjP ^2 

^2 — / =-^2;diag / 

^Jupui 

= Trel + Vff, ( 10 ) 

where is the resulting primary effective NN inter¬ 
action and we suppress the subscript “2.” The inter¬ 
action depends on A and the chosen P space in¬ 
cluding the selected value of ft. Note that the unitary 
transformation m is identical to OLS unitary transfor¬ 
mation which satisfies the decoupling condition 

QH^f^P := = 0 where the submatrix QH^f^P = 0 

is one of two decoupling conditions depicted in Fig. [1] for 
the primary Hamiltonian. 

There are certain freedoms within the OLS renormal¬ 
ization procedure as well as mathematical restrictions 
p^ . In this context, we note that in, our application, we 
select the dp lowest eigenvalues and eigenvectors of H 2 
for input to our primary effective Hamiltonian through 
Eq. © and obtain numerically stable and accurate re¬ 
sults. 


B. Transformation of many-body Hamiltonian into 
sd-shell space 


After a unitary transformation of the bare Hamilto¬ 
nian in Eq. (|4]) to the ATift (A^max = 4) model space 
for the case of ^^F, we calculate the 18-body effective 
Hamiltonian PH^^^P := in the ATift space and solve 
for its low-lying eigenvalues and eigenvectors in a NCSM 
calculation. This is analogous to solving the a = 2 case 
above so we introduce the corresponding subscript 18. 
We obtain a sufficient number of these 18-body solutions 
to generate a second unitary transformation to take 
from the ATift model space to a smaller secondary sub¬ 
space P', e.g., the sd-shell space, given by P^ax = 0- 
The secondary effective Hamiltonian is called P^ ^ with 
a' = 18 and is represented by P'P'^^P' in Fig. [Tj). 

This “second step” outlined above, follows a similar 
path to the “first step” and is indicated by the work flow 
at the upper right in Fig. [TJ Note that the P space in the 
first unitary transformation is now split into parts related 
to the two subspaces, P' and Q', where P' P Q' = P. 
Our secondary effective Hamiltonian Pf^^ is designed 
to reproduce exactly the lowest dp> eigenvalues of the 
primary effective Hamiltonian Pj^ through: 


P 


p'p 

18 


TjP t TjP' 

^18 ttP' ^18 

j ——--^ 18 ;diag /——- 

^jupup 

T,el + Vff ( 11 ) 


where ^ 


is the resulting secondary effective interac¬ 


tion and we suppress the label for the a' = IS depen¬ 
dence. 

This secondary effective Hamiltonian m is, in gen¬ 
eral, an 18-body operator. However, in the P^ax = 0 
case, the matrix dimension of the 18-body secondary ef¬ 
fective Hamiltonian m is the same as the matrix dimen¬ 
sion of a one-body plus two-body effective Hamiltonian 
acting in the sd-shell space. This means that Pf^ ^ can 
be taken to consist of only one-body and two-body terms, 
even after the exact 18-body cluster transformation. All 
the orbitals below the sd-shell space are fully occupied by 
the other 16 nucleon spectators, and the total 18-body 
wave function can be exactly factorized into a 16-body 
0+ and two-body sd-shell wave functions. This consid¬ 
erably simplifies calculations with Pfg^. Therefore, we 
can write a' as a' = Uc + av, where Uc is the number of 
core nucleons (16 in this case) and is the size of the 
valence cluster. 

In the third step outlined above we solve for the eigen¬ 
values of ^^F and in the P space using the effec¬ 
tive interaction from Eq. ([TQ|) joined with the Trei 

for A = 17, to obtain the proton and neutron one-body 
terms of the secondary effective Hamiltonian in the sd- 
shell space. Then we subtract the one-body terms from 
the secondary effective Hamiltonian of ^^E, and we ob¬ 
tain the effective “residual two-body interaction” matrix 
elements (or simply the TBMEs) in the sd-shell space. 
Additionally, in the third step, we evaluate the core 
energy by solving for its ground-state energy using the 
effective interaction from Eq. (pT|) joined with Trei 
for 1^0. 

Here, we adopt the A = 16 (17) relative-kinetic-energy 
operators for NCSM evaluations of the core (Angle- 
particle) energies in step three. In earlier papers [2l|, [^, 
based on the NCSM with a core first developed in 
Ref. a much stronger A dependence was obtained 
than in our present sd-shell calculations. We now under¬ 
stand these earlier results in terms of how the core and 
single-particle energies are calculated. In these earlier 
studies, the A dependence of the kinetic-energy opera¬ 
tor in the many-nucleon Hamiltonian used for calculat¬ 
ing the core and single-particle energies was taken to be 
the total A of the nucleus being studied. In our cur¬ 
rent calculations, we use A (core) = 16 for the kinetic- 
energy operator when calculating the core energy and 
A(core + 1) = 17 when calculating the single-particle en¬ 
ergies, independent of the total A of the nucleus being 
studied. Because the A dependence of the kinetic energy 
operator goes as 1/A, using the total A instead of A(core) 
or A (core -h 1) produces a much larger A dependence of 
the core and single-particle energies. Both these choices 
are technically correct (i.e., they produce identical re¬ 
sults for the nucleus being studied as we have verified), 
and merely reflect that these effective valence-space in¬ 
teractions are not uniquely defined. With our current 
choice, we achieve weak A dependence of our resulting 
core, single-particle and valence effective interactions for 
sd-shell applications, which is appealing since this is a 











5 


characteristic that is commonly found in phenomenologi¬ 
cal effective interactions. Our weak A dependence is also 
consistent with other ab initio investigations using either 
the IM-SRG technique [22| or the coupled cluster method 

We then proceed to the fourth step and subtract this 
core energy from the energies of the single-particle states 
of and mentioned above to arrive at our va¬ 
lence single-particle energies. At the completion of step 
four, we have our twobody valence-core (2BVC) effective 
Hamiltonian that may be used in standard shell-model 
(SSM) calculations. 

By using these core plus valence space single-particle 
energies along with the derived residual two-body effec¬ 
tive interactions, we can perform the SSM calculations for 
as well as other nuclei, in the sd shell and compare 
with full NCSM calculations in the Ahft space by using 
the primary effective Hamiltonian. The SSM calculations 
for ^^F will, by construction, give the same results as the 
NCSM calculations for ^^F within numerical precision. A 
corresponding approach for A > 18 nuclei is exemplified 
below where we also provide a direct comparison between 
NCSM and SSM results. One may then proceed, in prin¬ 
ciple, with SSM calculations to cases where full NCSM 
results are beyond current technical means. 

We may summarize the results of steps 2-4 by arrang¬ 
ing the results for the secondary effective Hamiltonian 
^ into separate terms: 


for the primary and secondary effective Hamiltonians. 
Therefore, the A^^ax ~ ^ secondary effective Hamiltonian 
not only reproduces the appropriate A^max NCSM eigen¬ 
values but also affords access to wave functions for these 
^max = 0 states which may be written with a factorized 
CM wave function of the entire system. 


III. EFFECTIVE TWO-BODY sd-SHELL 
INTERACTION 

In NCSM calculations, the dimension of the primary ef¬ 
fective Hamiltonian increases very rapidly as we increase 
A^max and/or the number of nucleons. We restricted the 
model space to Vmax = 4 in order to limit the com¬ 
putational effort, since our main goal is to demonstrate 
the procedure to obtain effective interactions in the sd 
shell for the shell model with the core using the ab 
initio NCSM and to test these derived effective interac¬ 
tions with SSM calculations. In order to carry out NCSM 
calculations, we used the MFDn code with the 

JISPI6 and chiral N3LO NN interactions. For the SSM 
calculations, we used a specialized version of the shell- 
model code ANTOINE 

A. Core and valence effective interactions for the 
A = 18 system 


uP'P _ uP'P 
11 ^/ — £ 1 ^ 




p'p 


sp 




p'p 


( 12 ) 


where we have allowed for the more general case of two 
successive renormalization steps (signified by P'P) with 
a' = A in the present discussion. In Eq. (H^ Ha^ rep¬ 
resents the core Hamiltonian for Oc nucleons; idgp rep¬ 
resents the valence nucleon single-particle Hamiltonian 
and Va^ represents the a^-body residual effective valence 
interaction. Note that rnay be used for systems with 
more than Oy valence nucleons, as we demonstrate be¬ 
low. We also note that the core and the valence single¬ 
particle Hamiltonians include their respective kinetic- 
energy terms. 

In line with our approximations mentioned above, we 
use ^^E alone to derive our isospin-dependent effective 
two-body interaction ^ for the sd shell. We then 
restrict our applications, at present, to cases with at most 
one proton in the sd shell. 

In SSM calculations, one typically uses only the i^sp 
and Va^ terms in Eq. m- In phenomenological Hamil¬ 
tonians Hsp is often taken from experiment and Oy = 2 
matrix elements are obtained by fits to properties of a set 
of nuclei. We will present detailed comparisons between 
our derived Hsp and terms with phenomenological 
interactions in a future presentation. 

There is an important distinction between our SSM 
calculations (with our Hamiltonian derived from the ab 
initio NCSM) and conventional SSM calculations with 
phenomenological interactions. We preserve the factor¬ 
ization of the CM motion throughout our derivation 


Eollowing the methods presented in Sec. H for H 2 
in Eq. m, we calculated the 18-body primary effective 
Hamiltonians with Vmax = 4 and hCl = 14 MeV by us¬ 
ing the bare JISPI6 and chiral N3LO [28| potentials 
for . We chose hft = 14 MeV since it is near the 
minimum of the ground-state energy of at Vmax = 4 
[23 and it represents a typical choice for derived effec¬ 
tive shell-model valence interactions (see, for example. 
Ref. [13 )• Euture efforts with primary effective Hamil¬ 
tonians derived in larger-Vmax spaces will be needed for 
meaningful analyses of the hft dependence of our results. 

We solve for the ^^E spectra in NCSM calculations 
with these primary effective Hamiltonians and present 
the lowest 28 eigenvalues in Table [H The corresponding 
NCSM eigenvectors for these 28 states in the Vmax = 
4 space are the eigenvectors dominated by Vmax = 0 
components. These 28 eigenstates correspond with the 
complete set of Vmax = 0 states in the sd shell. 

Eor each of these primary effective Hamiltonians H 2 
we then followed steps 2-4 above to calculate secondary 
effective Hamiltonians H^g^ as well as the resulting six 
valence single-particle energies ^ (three for neutrons 
and three for protons) and 63 valence two-body effective 
interaction matrix elements of ^ in the coupled JT 
representation. 

We now elaborate on the method of separating the sec¬ 
ondary effective Hamiltonian H^g^ into its components 
indicated in Eq. m- According to step 3 we first per¬ 
form separate NCSM calculations for ^^E and using 
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TABLE I: The NCSM energies (in MeV) of the lowest 28 
states Ji of calculated in 4/iQ model space by using 
JISP16 and chiral N3LO NN interactions with TiQ. — 14 MeV. 



T 

JISP16 


T 

N3LO 

V 

0 

-122.742 

V 

0 

-126.964 

3^ 

0 

-122.055 

3^ 

0 

-126.214 

0 ^ 

1 

-121.320 

0 ^ 

1 

-125.510 

V 

0 

-120.329 

5^ 

0 

-124.545 

2 ^ 

1 

-119.505 

2t 

1 

-123.974 

2 ^ 

0 

-119.011 

2t 

0 

-123.890 

V 

0 

-118.709 


0 

-123.077 

Oj 

1 

-118.410 

oj 

1 

-122.586 

2 ? 

1 

-117.211 

2 ? 

1 

-121.588 

3^ 

1 

-117.035 

4^ 

1 

-121.512 

4^ 

1 

-117.004 

Q + 

1 

-121.450 

3? 

0 

-116.765 

3t 

0 

-121.376 


0 

-113.565 

V 

0 

-119.658 


0 

-112.314 

4j 

0 

-118.656 

2 : 

0 

-111.899 

2 t 

0 

-117.950 

V 

0 

-110.357 

V 

0 

-116.106 

4^ 

1 

-109.625 

4? 

1 

-115.785 

2 ^ 

1 

-109.292 

2 ? 

1 

-115.407 

V 

1 

-108.752 

3+ 

0 

-115.309 

Q + 

O 4 

0 

-108.706 

V 

1 

-114.870 

2 ? 

0 

-108.485 

2 ? 

0 

-114.787 

le 

1 

-108.055 

V 

1 

-114.392 

2 ? 

1 

-108.041 

3? 

1 

-114.258 

Q + 

1 

-107.874 

2 | 

1 

-114.176 

36 

0 

-101.528 

3^ 

0 

-109.316 

V 

0 

-99.946 

V 

0 

-107.798 

0 ^ 

1 

-99.848 

2 ^ 

1 

-107.473 

2 ? 

1 

-99.607 

0 ? 

1 

-107.436 


the Hamiltonian consisting of the same from Eq. (pd)|) 
combined with Trei for A = 17. These two calculations 
provide total single-particle energies for the valence pro¬ 
tons and neutrons, respectively, that are expressed as 
matrix elements of ^ -b ^. 

tic 

We continue with the second part of step 3 to obtain 
the core energy (E'core) through a NCSM calculation for 
by using the Hamiltonian consisting of from 
Eq. m in combination with T^ei for The result¬ 
ing ground-state energy defines the contribution of 

^ to the matrix elements of ^ obtained in 

(Xc (Xc t>P 

the and calculations. The valence single-particle 
energies, the eigenvalues of are then defined as the 

total single-particle energies less the core energy. 

To obtain the TBMEs of the valence effective interac¬ 
tion V‘f we execute step 4 and subtract the contri¬ 
butions of the core and valence single-particle energies 
from the matrix elements of ^ to isolate V 2 ' ^ in 
Eq. (Ha). To be specific, we designate our valence single- 


TABLE II: Proton and neutron single-particle energies (in 
MeV) for JISP16 effective interaction obtained for the mass 
of A = 18 and A = 19. 



00 

T-1 

II 

A = 19 


Tcore = —115.529 

Tcore = —115.319 


15 3 

15 3 

ji 

2 2 2 

2 2 2 

e". 

-3.068 -2.270 6.262 

-3.044 -2.248 6.289 

4^ 

0.603 1.398 9.748 

0.627 1.419 9.774 


TABLE III: Proton and neutron single-particle energies (in 
MeV) for chiral N3LO effective interaction obtained for the 
mass of A = 18 and A = 19. 



Tcore 

A =18 

= -118.469 

Tcore 

= 19 

= -118.306 


1 

5 

3 

1 

5 

3 

Ji 

2 

2 

2 

2 

2 

2 


-3.638 

-3.042 

3.763 

-3.625 

-3.031 

3.770 


0.044 

0.690 

7.299 

0.057 

0.700 

7.307 


particle states by their angular momenta ji = |, |. 

Then, we define the contribution to the doubly reduced 
coupled-JT TBMEs (signified by the subscript JT on the 
TBME) arising from the core and one-body terms as 

= (Ecore + (f^) 

where ej represents the valence single-particle energy for 
the orbital with angular momentum j and the super¬ 
script, n (p), designates neutron (proton) for the energy 
associated with the (^^E) calculation, respectively. 

The resulting doubly reduced coupled-JT TBMEs of 
the valence effective interaction ^ are expressed as 

{jajbWV^f ^\\jcjd)jT 

= UajbWK'^ - HC - (14) 

By using the symmetries of the coupled-JT representa¬ 
tion, there are 63 unique TBMEs for which ja < jb- 

We confirm the accuracy of this subtraction procedure 
by demonstrating that SSM calculations with the derived 
core, one-body, and two-body terms of Eq. m in the sd- 
shell space reproduce the absolute energies of the lowest 
28 states of the 4:hft NCSM calculations for ^^E shown 
in Table [H 

The results for the core energy (Tcore) and valence 
single-particle energies (e^,ep for the JISP16 interac¬ 
tion are presented on the left-hand side of Table [H] for 
our leading example where the primary effective Hamil¬ 
tonian is derived for A = 18. The corresponding core en¬ 
ergy and valence single-particle energy results for the chi¬ 
ral N3LO interaction are presented on the left-hand side 
of Table [ml The valence single-particle energies clearly 
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reflect overall Coulomb energy shifts between NCSM cal¬ 
culations for and 

The resulting TBMEs of the secondary effective Hamil¬ 
tonian in Eq. m and of the valence effective in¬ 

teraction ^ in Eq. m are given in the seventh and 
eighth columns respectively of Tables [IVl and in the 
Appendix. The results of Table IIVI are obtained with 
the JISP16 NN interaction while those in Table [A] are 
obtained with the chiral N3LO NN interaction. 

These results for A = 18 with JISP16 presented in Ta¬ 
bles [n] and HV] (as well as the corresponding results with 
chiral N3LO in Tables [ml and El show the dominant 
contribution of Ecore to the diagonal TBMEs of the sec¬ 
ondary effective Hamiltonian as may be expected. 

When these Ecore contributions along with the one-body 
contributions are subtracted following Eq. m, the re¬ 
sulting diagonal matrix elements of Vf ^ fall in the range 
of conventional phenomenological valence nucleon effec¬ 
tive interactions. The nondiagonal TBMEs for A = 18 
shown in columns seven and eight of Tables and lAl 
remain unchanged by the subtraction process of Eq. dm 
as required by the Kronecker deltas in Eq. m- 

The resulting TBMEs of in column eight of Ta¬ 

bles [l^ and [A] (see tables in the Appendix) appear highly 
correlated, as shown in Eig. [2l indicating significant in¬ 
dependence of the valence nucleon interactions from the 
underlying realistic NN interaction. On the other hand, 
there is a noticeable dependence on the NN interac¬ 
tion seen in the spin-orbit splitting of the valence single¬ 
particle energies in Tables [IT] and im Eor both the split¬ 
ting of the ^ 5/2 and dg /2 orbitals, and the splitting of the 
Si /2 and the dg /2 orbitals, the JISP16 interaction pro¬ 
duces significantly larger results than the chiral N3LO 
interaction. This is most noticeable in the approximately 
30%, or 2 MeV, larger splittings of the d ^/2 and dg /2 or¬ 
bitals obtained with JISP16. 

Note that both JISP16 and N3LO lead to splittings 
of the d ^/2 and dg /2 orbitals that are larger than the 
phenomenological shell-model result which is based on 
experiment. In addition, the order of the calculated Si /2 
and ^ 5/2 orbitals are inverted compared with experiment. 
That is, ^^E and have a 5/2+ ground state, an ex¬ 
cited 1/2+ at about 0.5 and 0.9 MeV respectively, and an 
excited 3/2+ at about 5 MeV. Of course, neither JISP16 
nor the chiral N3LO interaction have been fit to any ob¬ 
servable in the sd shell. In addition, these calculated 
splittings should be sensitive to the 3NE, which is known 
to impact spin-orbit comling effects in p-shell NCSM in¬ 
vestigations jl. [TBI - IT^ . . 

B. Two-body valence cluster approximation for 
A = 19 system 

We now illustrate our approach for going to heavier 
nuclei by adopting the specific example of ^^E. In theory, 
we could proceed as with our application in the previous 
section, retain Uc = 16 and increase ay in pace with the 



FIG. 2: (Color online) Correlation between chiral N3LO and 
JISP16 TBMEs plotted in units of MeV. The 63 TBMEs are 
derived for A = 18 with the methods described in the text and 
are presented in the eighth columns of Tables IIVI and m Red 
circles (blue squares) represent T = 0 (1) matrix elements. 
The diagonal dashed line is the reference for equal matrix 
elements. The solid red line is a linear fit to the correlation 
points with the result y — 0.981x. The root-mean-square 
deviation between the two sets of TBMEs is 0.203 MeV. A 
plot of the A = 19 results in the tenth columns of Tables IIVI 
and lAl would be nearly indistinguishable from this plot. 


increase with A. Thus, for A = 19 we would derive ma¬ 
trix elements of an effective valence 3NE. However, this 
is not a practical path since there is no net gain over per¬ 
forming full NCSM calculations for each A en route to 
the secondary effective Hamiltonian. Instead, we present 
an alternative approximate path to heavier nuclei. 

Our procedure for going to heavier nuclei in the sd 
shell is to specify the sd-shell nucleus of interest with 
its value of A in the first step - the construction of the 
primary effective interaction of Eq. (pT|) . Then we 
define the two-body cluster Hamiltonian in Eq. (j4]) with 
this new value of A (A = 19 in our specific example) 
which is subsequently used to construct the primary ef¬ 
fective interaction. Next, we perform steps 2-4 as before 
with Oc = 16, Oy = 2, and neglecting effective many- 
valence-nucleon interactions: we perform ^^E, ^^E, 
and NCSM calculations with this primary effective 
interaction in order to extract the core energy, proton 
and neutron valence single-particle energies, and valence 
TBMEs. This is the 2BVC applied for general A. The 
generalization to = 3 (the 3BVC approximation) is 
straightforward but computationally demanding. Note 
that for A = 19 the 3BVC would correspond to a com¬ 
plete NCSM calculation. 

As an alternative, one may simply neglect any A de¬ 
pendence of the core energy, valence single-particle ener¬ 
gies, and valence TBMEs and perform SSM calculations 
throughout the sd shell with the effective shell-model in¬ 
teraction derived for ^^E. We also illustrate this choice 
below with the example of ^^E. 
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We now investigate the consequences of neglecting the 
induced 3NF and of neglecting the A dependence of V^. 
That is, we simply use the the derived core energy, va¬ 
lence single-particle energies, and valence TBMEs from 
the previous section in a SSM calculation of For com¬ 
parison, we also derive these quantities specifically for the 
^^F system in the 2BVC approximation, and we compare 
both with a complete NCSM calculations for ^^F, which 
corresponds to performing the 3BVC approximation. 

For the 2BVC approach to ^^F, we perform step 1 be¬ 
ginning with ^4 = 19 instead of ^4 = 18 in Eqs. 0-611. 
That is, we calculate the primary effective Hamiltonian of 
Eq. (pT|) for ^^F instead of ^^F. Then we proceed through 
the remaining equations, as we did for ^^F, using de¬ 
fined in Eq. (unj. For example, in the second step we 
solve for the secondary effective Hamiltonian ^ with 
a' = 18 at A^max = 4 using Eq. m as before. This es¬ 
tablishes the foundation for proceeding with steps 3 and 
4 to obtain the core energy, valence single-particle ener¬ 
gies and valence TBMEs needed for solving ^^F in a SSM 
calculation. 

The resulting core energies and valence single-particle 
energies calculated by using JISP16 and chiral N3L0 ef¬ 
fective interactions are given in the right-hand columns 
of Tables m and imi respectively. The core energies for 
the A = 19 case are less attractive than the ^4 = 18 case 
by 210 keV (163 keV) for JISP16 (chiral N3L0). The 
single-particle energies for the A = IS and ^4 = 19 cases 
differ by less than 30 keV (20 keV) for JISP16 (chiral 
N3LO). We observe, therefore, that the core and single¬ 
particle energies exhibit similarly weak A dependence for 
both interactions. 

The resulting TBMEs of the secondary effective Hamil¬ 
tonian in Eq. (pT]) and of the valence effective in¬ 

teraction V‘f ^ in Eq. (fT4|) are given in the ninth and 
tenth columns, respectively, of Table IIVI (for JISP16) 
and Table [X] (for chiral N3LO) in the Appendix. One 
observes a good correlation between the TBME results 
from the A = IS case and the A = 19 case by comparing 
column seven with column nine and column eight with 
column ten in both Tables [IVl (for JISP16) and (foi* 
chiral N3LO). The TBME’s of ^ exhibit particularly 
weak A dependence. The largest difference between the 
TBMEs in column eight and column ten in Table HVl ffor 
JISP16) is 9 keV and corresponding largest difference in 
Table CT Tor chiral N3LO) is 4 keV. 

Our observed weak A dependence of the core ener¬ 
gies, valence single-particle energies, and TBMEs is con¬ 
sistent with the view that the OLS transformation to the 
P space accounts for the high-momentum components of 
the NN interaction and the results are approximately 
independent of whether the two-body cluster is treated 
as embedded in A = 18 or in A = 19. The similarity of 
the derived TBMEs is also suggestive of a common, or 
universal, soft effective NN interaction. 

We may elaborate on these points by noting that the 
first OLS transformation can be viewed as reducing the 
ultraviolet (UV) regulator of the JISP16 and N3LO in¬ 


teractions to the UV scale of the HO basis space limit 
controlled by V^ax and by hCl. The HO basis space UV 
regulator imposed by our first OLS transformation may 
be estimated by using V, the maximum of 2n + / of the 
HO single-particle orbits included in the P space. For ^^F 
(or ^^F) with Vmax = 4 and hft = 14 MeV this UV reg¬ 
ulator is estimated to be either {N 3/2)mU = 1.59 
fm-i [al or ^2{N + 3/2 + 2)mn = 2.53 fm-i [H. In 
either case, the estimated UV regulator is independent of 
A and is sufficiently low that we may speculate that our 
chosen NN interactions are yielding a common (or uni¬ 
versal) UV-regulated primary effective W-interaction 
with the UV-regulation scale fixed by our choice of P 
space. Subsequent processing through the second OLS 
transformation is the same for both primary effective NN 
interactions so it retains that universality feature. 


C. SSM and NCSM calculations for and ^®F 
with A = 18 and A = 19 interactions 

We performed SSM calculations for the ground state 
and a few low-lying excited states of ^^F and ^^F by using 
the secondary effective Hamiltonians ^ of Eq. (H^ 
developed from the JISP16 and chiral N3LO potentials. 
We performed these SSM calculations with the code AN¬ 
TOINE by explicitly summing the one-body and 

two-body components on the right-hand side of Eq. m 
whose matrix elements are presented in Tables [IT] and 
llVlfor JISP16 and in Tables [ml and [A] for chiral N3LO. 
Next, we add the respective core energy to the result¬ 
ing spectra to yield total energies for comparison with 
NCSM calculations performed with the primary effective 
Hamiltonian. 

We also carried out NCSM calculations for ^^F and 
^^F by using the primary effective Hamiltonians H 2 of 
Eq. m which are based on the selected NN interac¬ 
tion and on the selected A in Eqs. ([I])-(|4]). The SSM 
and NCSM results for the ground state and a few low- 
lying excited states of ^^F and ^^F are shown in Fig. [3] 
for JISP16 and in Fig. 0] for chiral N3LO. The nucleus 
for which the spectra are presented (^^F or ^^F) is speci¬ 
fied at the top of each column along with the many-body 
method - either NCSM with the primary effective Hamil¬ 
tonian or SSM with the secondary effective Hamiltonian. 
Below each column we specify the A used in Eqs. 0 - 0 . 
When the results of the NCSM and SSM are the same 
with both many-body methods (as they should be the¬ 
oretically for ^^F), they appear as a single column with 
the label “NCSM/SSM.” This situation, a simple cross 
check of the manipulations and the codes, is presented 
in the first column of Fig. [3] for JISP16 and of Fig. lUfor 
chiral N3LO. Although these figures show only the lowest 
states, the cross-check is verified for all 28 states of two 
nucleons in the sd shell. 

The remaining three columns of Figs. [3] and 0] display 
two SSM calculations with the secondary effective Hamil¬ 
tonians and the exact NCSM calculation, all for ^^F. The 
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FIG. 3: (Color online) The ground-state energy (in MeV) and 
low-lying excited-state energies of and ^®F obtained by 
the NCSM and SSM calculations using the effective JISP16 
interaction. The tags A = IS and A = 19 at the bottom 
of each column refer to the effective JISP16 interaction ob¬ 
tained with the 2BVC approximation for general A. That is, 
the tags A = IS and A = 19 represent nucleus A used for 
deriving the primary effective Hamiltonian. In addition, we 
retain only effective core, one-body, and two-body terms for 
the secondary effective Hamiltonian. 
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FIG. 4: (Color online) The ground-state energy (in MeV) 
and low-lying excited-state energies of ^®F and ^®F obtained 
by the NCSM and SSM calculations using the effective chiral 
N3LO interaction. The tags A = IS and ^ = 19 at the bot¬ 
tom of each column refer to the effective JISP16 interaction 
obtained with the 2BVC approximation for general A. That 
is, the tags A = IS and ^ = 19 represent nucleus A used for 
deriving the primary effective Hamiltonian. In addition, we 
retain only effective core, one-body, and two-body terms for 
the secondary effective Hamiltonian. 
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second (third) column shows the results of using the pri¬ 
mary effective Hamiltonian for A = IS {A = 19) in the 
2BVC approximation and solving the resulting SSM for 
^^F as outlined above. The difference between the sec¬ 
ond and third columns is interesting since it reflects two 
different 2BVC approximations. In the second column, 
we see the effect of ignoring the contributions (both two 
body and three body) that one additional neutron makes 
by interacting with all nucleons in ^^F. In the third col¬ 
umn we see the effect of ignoring the contributions of all 
interactions in ^^F to the effective three-body valence in¬ 
teraction in ^^F. The differences between columns two, 
three, and four (discussed further below) are almost en¬ 
tirely due to the differences in the ground-state energies; 
the spectra are nearly the same. The ground-state en¬ 
ergies in columns two and three in Figs. [3] and m differ 
over a range from 4 to 211 keV compared with the exact 
results in column four. 

The effects neglected in the two different approxima¬ 
tions represented in columns two and three of Figs. [3] and 
mied to small differences in the spectroscopy and, there¬ 
fore, suggest that both are potentially fruitful paths for 
further investigation. However, when performing 2BVC 
calculations for M > 19 nuclei (i.e., continuing to retain 
only core-, one- and two-body interaction terms) it is nat¬ 
ural to expect that the difference between the SSM and 
NCSM calculations would increase due to the neglect of 


induced valence three-body, four-body, etc., interactions. 
The current results suggest that the dominant effect of 
neglecting these higher-body induced interactions may 
appear mainly as an overall shift in the spectrum. For 
the case of ^^F, the shift between columns three and four 
in Fig. [3] (j4]) shows that the 2BVC approximation for 
M = 19 is responsible for an overall net attraction (re¬ 
pulsion) of about 117 keV (4 keV) which is small on the 
scale of the overall binding. 

The overall shift between columns two and three in 
Fig. [3] dl]) shows that the differences in our derived SSM 
Hamiltonians produce about a 284 keV (207 keV) dis¬ 
placement in the binding energy. By referring to the 
results shown in Tables nn and imi we find that this 
displacement in binding energies is attributed approxi¬ 
mately to the difference in the core energies (about 80% 
of the displacement) and to the difference in the sum 
of single-particle energies for the three valence nucleons 
(about 20% of the displacement). These displacements 
may be cast either as diagonal matrix elements of ne¬ 
glected induced 3NFs or as corrections to the core and 
valence single-particle energies (or to a combination of 
both). The distribution of these displacements will ap¬ 
pear naturally when the 3BVC (i.e., full ay = 3) calcula¬ 
tion is performed for ^^F. 
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IV. SUMMARY, CONCLUSIONS, AND 
OUTLOOK 

We calculated A-dependent effective NCSM Hamilto¬ 
nians, called primary effective Hamiltonians [step 1 that 
leads to Eq. m], in a Ahft model space with the realis¬ 
tic JISP16 and chiral N3LO NN interactions. Next, we 
have solved the NCSM for low-lying eigenstates sufficient 
to derive a secondary effective Hamiltonian that acts only 
in the OhQ model space for the sd shell yet retains infor¬ 
mation from the full A-body correlations present in the 
NCSM solutions [step 2 that leads to Eq. (pA]) ]. We then 
separate the TBMEs of the secondary effective Hamil¬ 
tonians into core, one-body, and two-body contributions 
[steps 3 and 4 that lead to Eq. ([12])] which defines to 
the 2BVC effective Hamiltonian suitable for SSM calcula¬ 
tions. Einally, we use these secondary effective Hamilto¬ 
nians in SSM calculations and compare with exact NCSM 
results based on the primary Hamiltonians for A = 18 
and 19. 

We estimate that the first OLS transformation on the 
JISP16 and chiral N3LO NN interactions produces pri¬ 
mary effective interactions down to a sufficiently low 
UV regulator scale that we obtain a nearly common, 
or universal, primary effective NN interaction. Subse¬ 
quent processing through the second OLS transformation 
retains universality features resulting in TBMEs from 
JISP16 and chiral N3LO that are highly correlated as 
visualized in Eig. [2j 

The SSM spectra for A = 18 in the valence space are 
the same as the low-lying NCSM spectra since our the¬ 
ory of the secondary effective Hamiltonian is derived from 
the NCSM solutions obtained with the primary effective 
Hamiltonian. With the 2BVC approximation, for which 
we present two approaches, there are small differences 
between the SSM and NCSM spectra for the ^^E system. 
These differences are due to the omitted three-body ef¬ 
fective interactions for the ^^E system and are observed 
primarily as overall shifts in the spectra that are mainly 
due to shifts in the core energies. Close examination of 
the core, one-body, and two-body components of the sec¬ 
ondary effective Hamiltonians shows weak A dependence, 
which is encouraging for applications to heavier nuclei. 

We will extend our investigations to obtain more com¬ 
plete results in sd shell by proceeding to a higher V^ax 
model space for NCSM solutions with the primary effec¬ 
tive Hamiltonian. We will extend the 2BVC approxima¬ 
tion to the 3BVC approximation by including the three- 
body components of the secondary effective Hamiltoni¬ 
ans. In addition, we plan to incorporate initial 3NEs 
in the NCSM calculations that complement the realistic 
NN interactions. 
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Appendix A: Tabulation of derived two-body matrix 
elements 

In this Appendix, we present the tables of our derived 
2-body matrix elements (TBMEs). 

TABLE IV: The TBMEs (in MeV) of the secondary sd-shell 
effective Hamiltonian Hig ^ obtained from the NCSM calcu¬ 
lation with Amax = 4, TiVt = 14 MeV, and JISP16 potential 
for are shown as well as the TBMEs of its residual valence 
effective interaction, V‘f ^. Pairs of columns are labelled by 
the A used in Eqs. ffl-e to develop the primary effective 
NCSM Hamiltonian as discussed in the text. 


2ja 

2j6 

2jc 

24 d 

J 

T 

A = 

18 

A = 

19 

ttP'P 

^18 

vr'^ 

ttP'P 


1 

1 

1 

1 

0 

1 

-120.176 

-2.182 

-119.917 

-2.181 

1 

1 

3 

3 

0 

1 

-0.924 

-0.924 

-0.924 

-0.924 

1 

1 

5 

5 

0 

1 

-1.274 

-1.274 

-1.274 

-1.274 

3 

3 

3 

3 

0 

1 

-100.477 

-0.958 

-100.214 

-0.958 

3 

3 

5 

5 

0 

1 

-3.397 

-3.397 

-3.396 

-3.396 

5 

5 

5 

5 

0 

1 

-118.926 

-2.525 

-118.673 

-2.525 

1 

1 

1 

1 

1 

0 

-121.296 

-3.302 

-121.032 

-3.296 

1 

1 

1 

3 

1 

0 

-0.378 

-0.378 

-0.383 

-0.383 

1 

1 

3 

3 

1 

0 

0.231 

0.231 

0.236 

0.236 

1 

1 

3 

5 

1 

0 

2.054 

2.054 

2.052 

2.052 

1 

1 

5 

5 

1 

0 

-0.936 

-0.936 

-0.939 

-0.939 

1 

3 

1 

3 

1 

0 

-112.168 

-3.412 

-111.902 

-3.406 

1 

3 

3 

3 

1 

0 
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-1.380 
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-1.384 

1 

3 

3 

5 

1 

0 
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1.455 
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1 

3 

5 

5 

1 

0 
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3 

3 

3 

3 

1 

0 

-100.450 

-0.931 

-100.181 

-0.925 

3 

3 

3 

5 

1 

0 

-0.172 

-0.172 
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3 

3 

5 

5 

1 

0 
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3 

5 

3 

5 

1 

0 
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5 

5 

5 

1 

0 
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3.580 
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5 

5 

5 

5 

1 

0 

-117.448 

-1.047 
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-1.043 

1 

3 

1 

3 

1 

1 

-108.749 

0.007 

-108.487 

0.009 

1 

3 

3 

5 

1 

1 

0.042 

0.042 

0.042 

0.042 
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TABLE IV: {Continued.) 
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5 

3 
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-108.057 

-0.097 

-107.798 

-0.096 

1 

3 

1 

3 

2 

0 

-110.023 

-1.267 

-109.760 

-1.264 

1 

3 

1 

5 

2 

0 

-2.969 

-2.969 

-2.968 

-2.968 

1 

3 

3 

5 

2 

0 

-1.873 

-1.873 

-1.873 

-1.873 

1 

5 

1 

5 

2 

0 

-117.279 

-0.081 

-117.021 

-0.079 

1 

5 

3 

5 

2 

0 

-1.597 

-1.597 

-1.597 

-1.597 

3 

5 

3 

5 

2 

0 

-112.093 

-4.133 

-111.826 

-4.124 

1 

3 

1 

3 

2 

1 

-109.374 

-0.618 

-109.113 

-0.617 

1 

3 

1 

5 

2 

1 

1.504 

1.504 

1.504 

1.504 

1 

3 

3 

3 

2 

1 

0.185 

0.185 

0.185 

0.185 

1 

3 

3 

5 

2 

1 

0.601 

0.601 

0.601 

0.601 

1 

3 

5 

5 

2 

1 

1.005 

1.005 

1.005 

1.005 

1 

5 

1 

5 

2 

1 

-118.667 

-1.469 

-118.411 

-1.469 

1 

5 

3 

3 

2 

1 

-0.840 

-0.840 

-0.840 

-0.840 

1 

5 

3 

5 

2 

1 

-0.374 

-0.374 

-0.374 

-0.374 

1 

5 

5 

5 

2 

1 

-0.780 

-0.780 

-0.780 

-0.780 

3 

3 

3 

3 

2 

1 

-99.766 

-0.247 

-99.503 

-0.247 

3 

3 

3 

5 

2 

1 

-0.933 

-0.933 

-0.933 

-0.933 

3 

3 

5 

5 

2 

1 

-0.730 

-0.730 

-0.730 

-0.730 

3 

5 

3 

5 

2 

1 

-108.232 

-0.272 

-107.973 

-0.271 

3 

5 

5 

5 

2 

1 

-0.352 

-0.352 

-0.352 

-0.352 

5 

5 

5 

5 

2 

1 

-117.617 

-1.216 

-117.364 

-1.216 

1 

5 

1 

5 

3 

0 

-121.030 

-3.832 

-120.770 

-3.828 

1 

5 

3 

3 

3 

0 

0.068 

0.068 

0.066 

0.066 

1 

5 

3 

5 

3 

0 

1.373 

1.373 

1.375 

1.375 

1 

5 

5 

5 

3 

0 

-1.766 

-1.766 

-1.768 

-1.768 

3 

3 

3 

3 

3 

0 

-102.271 

-2.752 

-102.006 

-2.750 

3 

3 

3 

5 

3 

0 

2.000 

2.000 

1.998 

1.998 

3 

3 

5 

5 

3 

0 

0.961 

0.961 

0.963 

0.963 

3 

5 

3 

5 

3 

0 

-108.629 

-0.669 

-108.367 

-0.665 

3 

5 

5 

5 

3 

0 

2.308 

2.308 

2.306 

2.306 

5 

5 

5 

5 

3 

0 

-117.125 

-0.724 

-116.870 

-0.722 

1 

5 

1 

5 

3 

1 

-117.022 

0.176 

-116.765 

0.177 

1 

5 

3 

5 

3 

1 

-0.356 

-0.356 

-0.356 

-0.356 

3 

5 

3 

5 

3 

1 

-107.888 

0.072 

-107.629 

0.073 

3 

5 

3 

5 

4 

0 

-112.314 

-4.354 

-112.049 

-4.347 

3 

5 

3 

5 

4 

1 

-109.863 

-1.903 

-109.605 

-1.903 

3 

5 

5 

5 

4 

1 

-1.303 

-1.303 

-1.303 

-1.303 

5 

5 

5 

5 

4 

1 

-116.766 

-0.365 

-116.513 

-0.365 

5 

5 

5 

5 

5 

0 

-120.329 

-3.928 

-120.075 

-3.927 


TABLE V: The TBMEs (in MeV) of the secondary sd-shell 
effective Hamiltonian Hig ^ obtained from the NCSM cal¬ 
culation with iVmax = 4, %CL — 14 MeV, and chiral N3LO 
potential for ^^E are shown as well as the TBMEs of its resid¬ 
ual valence effective interaction, V-f ^. Pairs of columns are 
labelled by the A used in Eqs. 0-® to develop the pri¬ 
mary effective NCSM Hamiltonian as discussed in the text. 


2ja 


2 jc 

2 jd 

J T 


18 


19 

ttP'P 

^18 


ttP'P 

^18 

vr'p 

1 

1 

1 

1 

0 

1 

-124.196 

-2.106 

-123.978 

-2.104 

1 

1 

3 

3 

0 

1 

-0.991 

-0.991 

-0.991 

-0.991 

1 

1 

5 

5 

0 

1 

-1.268 

-1.268 

-1.268 

-1.268 

3 

3 

3 

3 

0 

1 

-108.265 

-0.858 

-108.086 

-0.857 

3 

3 

5 

5 

0 

1 

-3.538 

-3.538 

-3.537 

-3.537 

5 

5 

5 

5 

0 

1 

-123.099 

-2.278 

-122.914 

-2.277 

1 

1 

1 

1 

1 

0 

-125.152 

-3.089 

-124.960 

-3.086 

1 

1 

1 

3 

1 

0 

-0.022 

-0.022 

-0.023 

-0.023 

1 

1 

3 

3 

1 

0 

-0.175 

-0.175 

-0.175 

-0.175 

1 

1 

3 

5 

1 

0 

2.315 

2.315 

2.314 

2.314 

1 

1 

5 

5 

1 

0 

-0.750 

-0.750 

-0.750 

-0.750 

1 

3 

1 

3 

1 

0 

-118.632 

-3.870 

-118.418 

-3.866 

1 

3 

3 

3 

1 

0 

-1.149 

-1.149 

-1.148 

-1.148 

1 

3 

3 

5 

1 

0 

1.568 

1.568 

1.568 

1.568 

1 

3 

5 

5 

1 

0 

0.355 

0.355 

0.355 

0.355 

3 

3 

3 

3 

1 

0 

-108.280 

-0.873 

-108.101 

-0.872 

3 

3 

3 

5 

1 

0 

-0.217 

-0.217 

-0.217 

-0.217 

3 

3 

5 

5 

1 

0 

2.265 

2.265 

2.264 

2.264 

3 

5 

3 

5 

1 

0 

-119.761 

-5.620 

-119.549 

-5.616 

3 

5 

5 

5 

1 

0 

3.377 

3.377 

3.375 

3.375 

5 

5 

5 

5 

1 

0 

-121.832 

-1.011 

-121.646 

-1.009 

1 

3 

1 

3 

1 

1 

-114.838 

-0.076 

-114.627 

-0.075 

1 

3 

3 

5 

1 

1 

-0.157 

-0.157 

-0.156 

-0.156 

3 

5 

3 

5 

1 

1 

-114.478 

-0.337 

-114.268 

-0.335 

1 

3 

1 

3 

2 

0 

-116.128 

-1.576 

-116.126 

-1.574 

1 

3 

1 

5 

2 

0 

-2.623 

-2.623 

-2.622 

-2.622 

1 

3 

3 

5 

2 

0 

-1.980 

-1.980 

-1.978 

-1.978 

1 

5 

1 

5 

2 

0 

-121.972 

-0.503 

-121.757 

-0.501 

1 

5 

3 

5 

2 

0 

-1.703 

-1.703 

-1.702 

-1.702 

3 

5 

3 

5 

2 

0 

-118.482 

-4.341 

-118.271 

-4.338 

1 

3 

1 

3 

2 

1 

-115.422 

-0.660 

-115.211 

-0.659 

1 

3 

1 

5 

2 

1 

1.569 

1.569 

1.569 

1.569 

1 

3 

3 

3 

2 

1 

0.188 

0.188 

0.188 

0.188 

1 

3 

3 

5 

2 

1 

0.695 

0.695 

0.695 

0.695 

1 

3 

5 

5 

2 

1 

0.883 

0.883 

0.883 

0.883 

1 

5 

1 

5 

2 

1 

-122.903 

-1.434 

-122.688 

-1.432 

1 

5 

3 

3 

2 

1 

-0.869 

-0.869 

-0.869 

-0.869 

1 

5 

3 

5 

2 

1 

-0.298 

-0.298 

-0.298 

-0.298 

1 

5 

5 

5 

2 

1 

-0.802 

-0.802 

-0.802 

-0.802 

3 

3 

3 

3 

2 

1 

-107.666 

-0.259 

-107.487 

-0.258 

3 

3 

3 

5 

2 

1 

-0.885 

-0.885 

-0.885 

-0.885 
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TABLE V: {Continued.) 


2 / a 


Ud J T 


18 


19 

ttP'P 

^18 


ttP'P 

^18 


3 

3 

5 

5 

2 

1 

-0.813 

-0.813 

-0.813 

-0.813 

3 

5 

3 

5 

2 

1 

-114.549 

-0.408 

-114.340 

-0.407 

3 

5 

5 

5 

2 

1 

-0.359 

-0.359 

-0.359 

-0.359 

5 

5 

5 

5 

2 

1 

-122.077 

-1.256 

-121.892 

-1.255 

1 

5 

1 

5 

3 

0 

-125.266 

-3.797 

-125.049 

-3.793 

1 

5 

3 

3 

3 

0 

0.155 

0.155 

0.154 

0.154 

1 

5 

3 

5 

3 

0 

1.206 

1.206 

1.205 

1.205 

1 

5 

5 

5 

3 

0 

-1.648 

-1.648 

-1.647 

-1.647 

3 

3 

3 

3 

3 

0 

-110.003 

-2.596 

-109.822 

-2.593 

3 

3 

3 

5 

3 

0 

1.819 

1.819 

1.818 

1.818 

3 

3 

5 

5 

3 

0 

0.564 

0.564 

0.563 

0.563 


TABLE V: {Continued.) 


‘^ja 

2jc 

2 /d 

J 

T 


18 


19 

ttP'P 

^18 


ttP'P 

^18 


3 

5 

3 

5 

3 

0 

-115.233 

-1.092 

-115.023 

-1.090 

3 

5 

5 

5 

3 

0 

1.940 

1.940 

1.939 

1.939 

5 

5 

5 

5 

3 

0 

-121.768 

-0.947 

-121.583 

-0.946 

1 

5 

1 

5 

3 

1 

-121.476 

-0.007 

-121.262 

-0.006 

1 

5 

3 

5 

3 

1 

-0.094 

-0.094 

-0.094 

-0.094 

3 

5 

3 

5 

3 

1 

-114.287 

-0.146 

-114.078 

-0.145 

3 

5 

3 

5 

4 

0 

-118.684 

-4.543 

-118.472 

-4.539 

3 

5 

3 

5 

4 

1 

-116.134 

-1.993 

-115.924 

-1.991 

3 

5 

5 

5 

4 

1 

-1.319 

-1.319 

-1.318 

-1.318 

5 

5 

5 

5 

4 

1 

-121.190 

-0.369 

-121.006 

-0.369 

5 

5 

5 

5 

5 

0 

-124.545 

-3.724 

-124.358 

-3.721 
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